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Abstract 

We present a semiclassical trace formula for the canonical partition function of arbitrary one- 
dimensional systems. The approximation is obtained via the stationary exponent method applied 
to the phase-space integration of the density operator in the coherent state representation. The 
formalism is valid in the low temperature limit, presenting accurate results in this regime. As 
illustrations we consider a quartic Hamiltonian that cannot be split into kinetic and potential 
parts, and a system with two local minima. Applications to spin systems are also presented. 
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I. INTRODUCTION 



The quantum partition function Z is the basic physical quantity from which all ensemble 
averages are derived in the canonical formalism p], 3, [J. It is given by the trace of the 
density operator Z = Tr [p], that can be expressed either as a discrete sum in the microstate 
energies Z = '^2 n e~ l3En or as a continuum integral, Z = 1/h J dqdp p, where p is some 
representation of the density operator in phase-space, like the Wigner function ^| or the 
coherent state average of p. In classical mechanics the density is simply p c i ass = e~^ n , where 
7i is the classical Hamiltonian. 

Because of the difficulties in the exact evaluation of Z, a variety of formal developments 
have appeared in the literature to include quantum effects in statistical mechanics, without 
exactly solving the Schrodinger equation. Semiclassical methods, in particular, have at- 
tracted a lot o attention P, UQ, 0, 0, Q • From a more basic point of view, it is well-known 
that classical statistical mechanics is, in most of cases, unable to reproduce the experimental 
behavior for T — > 0. A paradigmatic example is the failure at low temperatures of Dulong- 
Petit law for the specific heat of solids, or more generally, the failure of classical systems in 
obeying the third law of thermodynamics HQS . Roughly speaking, this is due to quantum 
discreteness that becomes relevant in this regime. Given this strong discrepancy between 
classical and quantum statistical mechanics for T — > 0, it is reasonable to ask what happens 
in a semiclassical formulation, since it should be somewhere in between the two theories. 
Therefore, it is desirable to obtain an explicit semiclassical expression for the canonical par- 
tition function, from which thermodynamic quantities such as internal energy and entropy 
can be derived in a closed analytical form in terms of classical ingredients. 

A straight, but not particularly effective, procedure in this direction is to replace the 
exact quantum levels by semiclassical ones in the discrete summation. For 1-D systems, 
expressions for the energy levels in the WKB approximation are relatively easy to obtain, 
but the corresponding partition function is inaccurate at low temperatures. More elaborated 
examples are the Wigner-Kirkwood method, which consists of an expansion of the density 
matrix in powers of Hj3 = H/IcbT approximate path integral representations and 
phase-space sampling with semiclassical quasi-probability distributions 0,0,0]. All these 
methods use classical ingredients to represent the density operator, but the procedure pre- 
sented in |f| is perhaps the most direct, in the sense that it simply replaces the classical 
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Boltzmann weight by a semiclassical one. In this paper we shall consider a similar method. 
The procedure is based on the recognition that the evolution operator in quantum mechanics 
and the Boltzmann operator in statistical physics are formally similar. Explicitly we have 
the correspondence e~ lHt ^ h — > e _/3W for t — > — it, with r = Hj3 (the so-called Wick rota- 
tion). The density operator p is obtained when such a rotation is applied to the diagonal 
quantum propagator and the partition function is written as a phase-space (or configuration 
space 0,13]) integral that can be calculated numerically. The ensemble averages are, con- 
sequently, also expressed in the form of integrations (see equations (2.16) to (2.18) in Q|). 
This method, in spite of its usefulness in numerical calculations, does not allow for a de- 
tailed analysis of the functional dependence of thermodynamic potentials, and other related 
physical quantities of interest, with the temperature. In this work we show that, in the low 
temperature limit, the phase space integrals appearing in the semiclassical expressions can 
actually be performed analytically. The resulting representation of the partition function in 
terms of classical quantities is similar to a trace formula, involving only equilibrium points 
and periodic orbits. In order to obtain the semiclassical partition function we shall use 
the coherent state representation. This choice enables the formalism to be applied in more 
general situations involving internal degrees of freedom such as spins. 

We believe that such a direct semiclassical formula should be of use in theoretical studies 
and also in numerical applications. Besides, some traditional methods like the Wigner- 
Kirkwood expansion converge very slowly for small T, justifying the particular interest in 
this region of temperatures. As we shall see, however, the semiclassical calculation of Tr[p] 
naturally leads to the complicated problem of summing over all periodic orbits of a related 
classical system. This is a common difficulty in practical applications of semiclassical trace 
formulas. Fortunatelly, due to the effects of the Wick rotation in the classical equations of 
motion, we shall see that this problem is not present in some important situations, e. g., 
one-dimensional potentials with a single minimum. 

The paper is organized as follows: in section II we review the concepts of Wick rotation, 
coherent states and semiclassical propagator. Next we derive the semiclassical formula for 
the partition function and give explicit expressions for potentials with a single well and 
multiple wells. In section IV we apply these formulas for three specific systems and compare 
our results with exact quantum and classical calculations and with other approximations. 
In section V we outline the extension of the formalism to spin systems. Finally, section VI 
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is devoted to some concluding remarks. 



II. PRELIMINARY DEFINITIONS 



The partition function can be written in terms of the quantum density of particles in the 
coherent state representation as 

Our procedure to calculate the semiclassical limit of Z consists of two main steps: first 
we perform a Wick rotation in the quantum propagator (z \ e~ %m / h \ z). Next we replace the 
propagator by its semiclassical expression and calculate the integral over q and p explicitly 
by the saddle point method to obtain Z sc . In the following subsections we review the basic 
ingredients of this approach: the Wick rotation, the coherent state representation and the 
semiclassical propagator. In the next section we calculate Z sc . 



A. Wick rotation and classical dynamics 



Semiclassical formulas contain classical quantities, such as actions and their derivatives, 
thus our first step is to understand the effect of the rotation t — > —it on these classical 
ingredients. We start with Newton's second law. Consider the pair of transformations 



t = -iT, V{q) = -V{q) 



(2) 



tfq 
dt 2 



(3) 



leading to 

dV .._dV _ .._ _dV 
dq ^ dq ^ dq 

where dot means derivative with respect to r. We see that in order to get back the usual 
form of the equation after the Wick rotation we have set V = —V, which means that the 
classical dynamics with complex time is formally equivalent to real time dynamics with the 
reversed potential V . The effect on the classical action is as follows (m = 1) 



S 



dt'C 



dt' 



-v 

2 1 dt 



ftp /i 



V 



lid 



dr'C = tS, (4) 



where it is understood that L and S in the right-hand side are evaluated for a final time 
hp and with the reversed potential V. There is, however, a more general approach, that we 



may call 'canonical Wick rotation', which is constructed in phase-space and does not rely 
on the existence of the function V(q). Starting from Hamilton's equations and making 



-IT 



H — iH (5) 



we obtain 



dq_dH dp _ dU _^ q --i dU p-i— _> q - dH p-- dH (6) 
dt dp ' dt dq ' dp ' dq dp ' dq 

which preserve the original form of the canonical equations. If the Hamiltonian can be 
written in 'Euclidean' form (p 2 /2 + V) the operations (j2J) and (jSJ) produce the same overall 
result. However, for more general Hamiltonians, for example with force fields not derivable 
from a potential or involving spin degrees of freedom, only the last pair of transformations 
is applicable. The classical action is invariant under these transformations, 
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In deriving our semiclassical partition function we shall employ the canonical Wick rotation 
because of its generality. As illustrations we apply the transformation (JHJ) to a non-Euclidean 
Hamiltonian and the transformation (0) to a system with a well defined potential in section 
IV. 



B. Coherent state representation 

Here we briefly describe some properties of the coherent state (or Bargman) representation 
of quantum mechanics. Canonical coherent states {\z}} can be expressed as an infinite sum 
of harmonic oscillator number states {|^)}: 

\ Z ) = e-'*l 2 / 2 y 41 |n) , where z = ±= (\ + i . (8) 

The real numbers q and p are the expectation values of the corresponding quantum operators 
q and p in the state \z), while the parameter b is proportional to the uncertainty in position, 
b = \f2Aq, and is related to the frequency of the associated harmonic potential by b 2 = 
h/muj. The resolution of unity can be written in terms of the set {\z)} as 
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In the coherent state representation the density of particles in phase-space reads 

p(q,p,(3) = (z\e-^\z) , (10) 

whose corresponding quantity under the inverse transformation r = it is the diagonal form 
of the coherent state propagator K(q,p,t) = (z\e~ tHt ^ n \z). 

We finish this subsection by defining the canonically conjugated variables u(t) = [q(t)/b+ 
i bp(t)/h]/ \[2 and v{t) = [q(t) /b — ibp(t) /h]/y/2 in terms of which Hamilton's equations read 

dv i &H du i &H 

dt = h~du ' dt = ~h~dv ' ( ' 

These variables are particularly convenient to write down the semiclassical limit of the 
coherent state propagator. 



C. Semiclassical coherent state propagator 

The set of coherent states forms a non-orthogonal over-complete basis, since each state 
in the set can be written as a linear combination of the others. This over-completeness 
implies the existence of several forms of the path integral formulation for the propagator, all 
equivalent quantum mechanically, but each leading to a slightly different semiclassical limit. 



Klauder and Skagerstam 12j proposed two basic forms for the coherent state path integral. 

The semiclassical limit of these propagators were considered in (kJ where it was shown that 

both propagators can be written in terms of classical complex trajectories, each governed 

by a different classical representation of the Hamiltonian operator 7i: the P representation 

in one case and the Q representation in other. The phase appearing in these semiclassical 

formulas turns is not just the action of the corresponding complex classical trajectory, but it 

also contains a 'correction term' that comes with different signs in each formula. In [13| it was 

also suggested that a semiclassical representation involving directly the Weyl representation 

of 7i, or the classical Hamiltonian Tl, could probably be constructed, and a formula for this 

representation was conjectured. This conjecture, along with the corresponding quantum 

mechanical path integral representation, has been recently proved ^| using the translation 

I I 

and reflection operators studied in In this paper we shall adopt this later semiclassical 
expression, since it is the simplest (although not always the most accurate |3]) of the three 
known formulas. This semiclassical expression for the propagator connecting an initial state 



I z') to a final state \z") is given by 



n\-^Ht\j 



[z e i 



(12) 



where there is an implicit sum over the, usually complex, trajectories satisfying Hamilton's 
equations (jll)) with boundary conditions u' = -u(O) = z', v" = t>(t) = z"*. The Hamiltonian 
entering in formula (j!2)l is the function 7i describing the corresponding classical system, and 
the function appearing in the exponent is the complex action given by 



o L 



ih ( du dv\ 

— v- u— — n 

2 \ dt dt J 



ih 



dt' - y [«V + V" U " 



(13) 



where we have set u" = u(t) and f' = f (0). Note that, due to the complex character of the 
classical orbits, in general, we have v' ^ z'* and u" ^ z". Finally, the pre-factor in formula 
(fT2"|) is written in terms of the tangent matrix, whose elements are given by the following 
relation 



\ 5v "; 



\ 



M vu M vv 



\ 5V '/ 



(14) 



where (Su', 5v') denote small initial deviations from the classical orbit and (Su", Sv") the 
corresponding deviations after a propagation time t. It is simple to show that det[M] = 1 
and that the following relations are valid 



ih 



-d 2 <b/du'dv'' 



d 2 §/du'dv" i 
d 2 §/dv" 2 ' M, 



d 2 $/du'dv" 
d 2 §/du' 2 



(15) 



III. THE SEMICLASSICAL CANONICAL PARTITION FUNCTION 



We proceed to make the transformations t = —ir and 7i = iH in Eq. fl2|) to get the 
semiclassical density of particles in phase-space. Since the complex action and, consequently, 
the tangent matrix elements remain invariant, the semiclassical density of particles at a given 
temperature is 



Ps 



[z e 



-m\j 



-i- expU$KX,T)-^V 
M vv I h 



(16) 



We remark that Hamilton's equations (|TT|) become 



with the boundary conditions (for the diagonal propagator) 



u (0) = u ' = z' , v{t) = v" = v{Hf3) = z' : 



(18) 



The semiclassical partition function is given by 

Z sr ~ 



dv! dv" 1 [ i , .. . . ... 

^ryw:, exp U^ v > u ^"=*^=*' ~ zz > ■ 



(19) 



In order to evaluate this integral by stationary exponent approximation we shall assume the 
low temperature limit (3 » 1. The stationarity conditions are 



f /* r\ f ft I n 

— Z — V — Z =0, — — — ; Z = U — Z = 



Hdu' 



H dv" 



where we used the relations d^/du' = —ihv' and d$>/dv" = —ihu" |13[. We get 



v' = z'*, u" = z'. 



(20) 



(21) 



Therefore, the trajectories must satisfy the four conditions (118)1 and (|21|). It can be shown 
that these restrictions demand the contributing orbits to be real and periodic. Next we 
expand the exponent in ()19)) around the stationary trajectories, labeled with the subscript 
'0', up to second order and evaluate the pre-factor on the contributing solutions. We obtain 



L <j>_„ V 

e« 



d(Au') d(Av") 

27H 



(22) 



1) Au'Av" 





where Au' = u' — u' and Av" = v" — v'q. The Gaussian integration gives 

-1/2 



v 7 ^ 



z <9 2 $ 
ft du'dv' 1 



+ 



1 <9 2 $ <9 2 $ 



^2 ^ M /2 ^ 



//2 



(23) 



From relations ()T5j) we obtain the following compact expression for the semiclassical partition 
function 



E 



g fi ^0 " " 



(24) 



, A /Tr[M,; j) ]-2 

where the sum runs over the real periodic trajectories with period hj3 and Tr[M] = M uu + 
M vv . Note that, since the contributing trajectories are periodic, the second term in the right- 
hand side of expression (JT3j) is exactly —ihu'v", so we can write the exponent in equation 
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as 

i 

n 



where we have used the periodicity in the last equality. Finally, we can write the expression 
for the semiclassical partition function as 



Z sc = Tr[p sc ] = £ Z { £ = £ / = • ( 26 ) 



We remark that, as is usual in stationary phase evaluations, there may exist divergent 
contributions corresponding to spurious stationary orbits. These orbits must not be taken 
into account in the semiclassical formula since their inclusion through the deformation of the 
original integration contours in (|19|) could not be justified by Cauchy's Integral Theorem. 
Note that the expected divergence in the infinite temperature limit occurs for r — ► 0, when 
M becomes the identity matrix and Tt[Mq ] — ► 2. Finally, we recall that, had we used 
transformation (J2J) instead of (|HJ) the above formula would have —S/h as the argument of 
the exponential. This alternative form is used in the next subsections. 



A. Single Well Potentials 

Now let us assume that there is a potential function V which has a single minimum 
(corresponding to a maximum of V). Due to the inversion in concavity, there are no periodic 
trajectories, except for the trivial orbit that, without loss of generality, we place at u = v = 
(or q = p = 0). So, in this important case, the summation in (j2E| reduces to a single term 
calculated on the trivial orbit: 

Z sc = , . (27) 

^TrfMo] - 2 

Formula (J2T)) can be evaluated more explicitly. We have the following relations for the time 
evolution of small perturbations: 

d 2 V 

5q = 5p, 5p=-—5q. (28) 

For a non-trivial trajectory the coefficient d 2 V/dq 2 = V' is a time dependent function, since 
it must be evaluated on the trajectory q(r), but for the trivial orbit it is simply a constant 
calculated at the maximum of V, and the above equations can be readily integrated. The 



result is Tr[M ] = 2 cos( V V r) = 2 cosh(VK" r), where V" > 0. This leads to the following 
semiclassical partition function 

e -(3V(0) 

2 SC = -. ;— . (29) 

2sinh UhVV"p) 

The above expression is simply the partition function for the harmonic oscillator that best 
fits V at q = 0. This result turns out to be too poor and the conclusion is that a higher 
order approximation is needed in this case. In section IV we illustrate such a higher order 
expansion for a specific example. In contrast, we shall show that the second order semiclassi- 
cal partition function given by (|2fij) already gives non-trivial results in the case of potentials 
with multiple minima. 



B. Multiple Well Potentials 

Let us assume that the potential V has N minima, located at q a , with a = 1, 2, . . . , N, 
and N — 1 local maxima. In this case there are non-trivial periodic orbits between two 
successive maxima of V. In addition each equilibrium point corresponds to a trivial orbit 
that, in principle, should contribute to the evaluation of Z sc . However, the N — 1 stable 
equilibrium points of V give rise to terms proportional to 1/ sin f^Vv^ ' (3\ , V' > 0, which 
have an oscillatory and divergent behavior for (3 — > oo, and correspond to the spurious 
saddle points we have mentioned. Discarding these contributions we get 

Z sc = V t r- + V . , (30) 

^^sinh^^v^t/?) " ^TrfM^l-2 

where we have split the sum in the trivial unstable equilibria (with respect to V) and the 
non-trivial periodic orbits (denoted by the primed sum). We call the first term the harmonic 
contribution and the second term the tunneling contribution, since the /z-orbits connect suc- 
cessive minima in the original potential. The calculation of the tunneling terms requires a 
careful procedure. For a fixed temperature T we must seek, in each well of V, orbits with 
period r = h/ksT and sum their contributions. Note also that these orbits should, in prin- 
ciple, contribute to the partition function at temperatures T/2, T/3,..., T/n, corresponding 
to the propagation times 2r, 3r,..., nr. However, these multiple traversal period orbits have 
larger and larger actions for increasing values of n, causing their contribution to fall off 
exponentially. 
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C. Connection with thermodynamics 



The connection between the canonical ensemble in statistical mechanics and thermody- 
namics can be made through relation Z = exp(— f3f), where / is the Helmholtz free energy 
per particle, whose semiclassical expression reads f sc = —kBTln(Z sc ). The internal energy 
can also be easily expressed in terms of classical quantities. Using the general formula (|26|) 
we obtain 

- = "I ,nZ - = h £ z « K + m HTllM ° )] ~ 4 ' (31) 

where we have used dS/ dr = —H = iH. Once we get the Helmholtz and internal energies we 
can determine the semiclassical expression for the entropy through relation s = (u — f)/T = 
&b/3(w — /). In the next section we give an example of a system that classically violates 
the third law, and show that, on the contrary, the semiclassical entropy does satisfy the 
condition s — > for T — > 0. Higher order derivatives of thermodynamical interest can also 
be obtained analytically, e. g., the specific heat du sc /dT: 



n ° ] + \§p HTT[M ° )] " 2) } 2 " lw HTr[M ° U)] 



(32) 

In the next section we shall apply our formalism to simple examples and compare the results 
with the pure classical and quantum ones. 

IV. APPLICATIONS 

A. The simple harmonic Hamiltonian 

We first illustrate our semiclassical procedure with the simple harmonic oscillator. The 
Hamiltonian is 

7i = -p 2 + 2 U}2< ^ 2 = fi JJJUV H — —ifiujuv . (33) 

The solutions of the equations of motion (|17jl are u(t') = -u'e - ^', v(r') = v " e a; ^ T '~ T - ) , with 
u(0) = z' and v(t) = z'*. From these relations we can write u" = u' e~ UJT and v" = v' e^, or 
5u" = e'^du' and 5v" = e UT 5v' . Therefore 

M uu = e~" T M m = e" T . (34) 
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The only periodic trajectory is u — v — 0, for which the complex action vanish. Gathering 
the ingredients together we get the exact quantum result for the partition function 

2^ \ (35") 

sc ~ Ve^P + e-^P - 2 ~ 2 sinh ' 1 1 

This expression leads to the well-known results for an ensemble of non-interacting harmonic 

oscillators. In particular, we get for c sc the Einstein expression for the specific heat of a 

crystalline solid (c sc = c einst ein)- 



B. The quartic non-Euclidean Hamiltonian: quadratic approximation 

We now consider a more challenging system described by the Hamiltonian 

n = ~L {l p2 + \ u2q2 + a ) 2 = m«« + ^ ' w 

with a = X/fkv. With the Wick rotation we get 

H = —ifhj(uv + a) 2 . (37) 

Note that this Hamiltonian may be seen as a sum of a harmonic term (proportional to a) 
and an anharmonic quartic term (independent of a). Since a ~ fT 1 , while uv ~ H°, in 
the semiclassical limit (a » 1) the quartic term must be seen as a perturbation added to 
the harmonic Hamiltonian. Note also that the H can not be written as p 2 /2 + V(q). The 
equations of motion are 

v — 2cu(uv + a)v , and u = — 2u(uv + a)u . (38) 

Since (uv + a) — ^/H/fko is a constant of the motion, the classical solution satisfying the 
boundary conditions is 

u(t') = u' e- nr ', v(t') = v" e^ T '' T ^ =}► u" = u' e~ n \ v" = v' e Ur , (39) 

where Q = 2u(u'v' + a) and r = hf3. We note that the only periodic orbit is the trivial one, 
although the system does not involve a potential. The above solution is very similar to that 
of the harmonic Hamiltonian but, in the present case, the frequency Q is energy dependent. 
The connection between initial and final variations in a given trajectory can be obtained 
from 

Su" = e- nT 5u' - Te' QT 5Q , 5v" = e QT 5v' + re nT 5Q , (40) 
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where 5CI = 2to(u'Sv' + v'Su'). Therefore the relevant tangent matrix elements are 

M uu = (1 - 2ujtu'v') e~ nT , (41) 

and 

M vv = (1 + 2ujtu'v') e nT . (42) 
On the trivial trajectory u = v = we obtain simply 

M uu = e" 2 ^, M vv = e 2a ^, ~S Q = -a 2 hw(5 , (43) 

where we have set Q = 2ua. This leads to the semiclassical result 

e -a 2 tkuf3 

Zsc = 2sinh(a^/?) ' (44) 

from which one can write 

f sc = a 2 ftw + — In [2 smh(afku {3)] , (45) 

and 

u sc = o^fiw + fkoa coth(aHujP) . (46) 

Not surprisingly, these expressions correspond to the exact quantum result for the harmonic 
part of the Hamiltonian (jHfijl . In the next subsection we develop a more accurate approxi- 
mation. The exact quantum partition function is given by Z = ^2 n e~^ En , where {E n } are 
the (non-degenerate) eigenvalues of the quantum Hamiltonian Tt = fijjj(a)a + a) 2 . Thus 

2 = e -^/ 3 (»+«+ 1 /2) 2 _ (47) 

n 

The classical function is given by 
where Erf denotes the error function. 



C. The quartic non-Euclidean Hamiltonian: higher order approximation 

Higher order corrections to the semiclassical partition function can also be obtained from 
(Jinj). We illustrate the method for the Hamiltonian (jHEJ)- Let us rewrite integral (fTT3|) as 

Z sc = j dU 2 ^ A exp u', t) v » =2 >* >u >= z > -z'vj , (49) 
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where A = 1/ y/M vv . Our procedure now is to expand the argument of the exponential up 
to fourth order around the saddle points, and the pre-factor up to second order (for more 

n 

details see appendix B of |13l|). We get 

A » A + A {u) Au' + A {v) Av" + -A iuu) Au' 2 + -A (vv) Av" 2 + A {uv) Au'Av" , (50) 

and 

i$ - z'*z' « - v? v» + f n <5> { uu)Au' 2 + ^ {vv) Av" 2 + (|$ M - l) Au'Av" 
+ -k [$ { uuu)Au' 3 + 3 <$> {uuv) Au' 2 Av" + 3 <5> {uvv) Au'Av" 2 + $ (wb) Aw" 3 ] 
+ 2ift [*(««u)A?i' 4 + 4 $ ( „) Au' 3 Ai;" + 6 $ (uuw) Am' 2 At;" 2 
+4 $( mm )Ait'At;" 3 + $ (ww) Aw" 4 ] , 

where A( u ) = (dA/du')o and so on. For the Hamiltonian ()36|) we find that the only non- 
vanishing terms are the ones proportional to A( uv ), and &( U uw) besides the zeroth 
order terms (see appendix). Therefore the partition function integral reads 

= e ^ j d(An0d(A,») (Aq + A{uv)Au > Av ») x 

exp { (i$ M - 1) Au'Av" + f h $ {uuvv) Au' 2 Av" 2 } 

(51) 

* e— 2h " J d ^ Av '' ] ( A o + A (uv) Au'Av") x 

(1 + ^ ®(uuw)Au' 2 Av" 2 ) exp { (i$ (wj) - l) An' Aw"} , 

where we have assumed that the fourth order contribution in the exponent is small if com- 
pared to the second order term, and the tilde indicates the higher order approximation. The 
next step is to diagonalize the quadratic form in the argument of the exponential. This can 
be done with the variables W and K given by 

Au' = -j={K-iW) , Av" = --—(K + iW) , (52) 

where 27 = {i^( uv \/h — 1). Integral (JoTj) becomes 



% 7 J 2ixi 



Ao - ^{W 2 + K 2 ) + (W 2 + K 2 ) 2 

27 lb nry z 



e -w 2 -K 2 



(53) 

where we have excluded terms of order higher than four. The resulting partition function is 
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where Z sc is given by expression (jUj). We have find that (see appendix) 

A = e~ a ^, A {uv) = -2hwf3e-* a ^ , <$> {uuvv) = Aih 2 uf3 e~ Aa ^ , (55) 
and, therefore, 

^"dSwC 1 -^^^ 1 ")- (56) 

In figure^we compare the scaled thermodynamic potentials /* = f/fvu and u* = u/Hu as 
functions of the reduced temperature T* = ksT/huj, with a = 8.0, for the four cases: higher 
order semiclassical (black), quantum (gray), classical (crosses), and quadratic semiclassical 
(dashed line). As expected the classical results tend to the exact quantum ones for high 
temperatures, while they completely fail to describe the physical behavior at T — > 0. In 
contrast, the semiclassical approximations are very precise in the low temperature limit. In 
particular we note that the semiclassical entropy satisfies the third law of thermodynamics 
while the classical function does not. This is shown in figure El where the scaled entropy 
s* = s/ks is plotted function of T*. 



D. Two-Well Potential 

Here we address a system with a quartic potential function V = A q A — B q 2 + C, with 
A, B > 0. If we denote the height of the local maximum located at q = by AE and the 
distance between the two symmetric minima by 2a, the potential is written as 



V(q) = AE 



q -Y-2( q -)\l 



(57) 



-a/ \aJ 

such that V(±a) = 0. In this case we apply equation (fHU|) directly. The harmonic contribu- 



tion is 2 x 1/ [2 sinh(ft\/ 2AEf3 / a)] , while the tunneling term must be determined numerically. 
There are two distinct regimes to be considered. If AE is large enough to encompass several 
energy levels, the harmonic contribution should be dominant over the tunneling one, since 
the two minima hardly 'see' each other. In this case the trajectories connecting the wells 
have large actions, giving exponentially small contributions to the semiclassical partition 
function. In contrast, when AE is comparable with the energy of the ground state, the 
tunneling term should give a relevant contribution. These statements are verified in the 
numerical computation of Z sc for AE = 3.0 and AE = 0.15 with a = 5 in both cases (we 
use arbitrary units for which H = 1 and ks — 1)- In the first case (deep well) the ground 
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state has an approximate energy of 0.5, and there are eight levels (four doublets) with en- 
ergy smaller than AE, as shown in fig. 0(a). Figure 0(b) shows the quantum partition 
function and Z sc , where we have used only the harmonic term, since the tunneling contri- 



bution is about three orders o 
and Patrascioiu (see also 



magnitude smaller. The approximation derived by Gildener 
) specifically for the potential (JH7|) is also displayed. This 
approximation and the semiclassical one produce almost identical results in this case. In 
the case of the shallow well there are only two levels with energy below AE = 0.15 (see 
fig E(a)). As expected the harmonic contribution is unable to give a good approximation 
to Z and the tunneling term must be considered. The overall semiclassical result is in very 
good agreement with the exact one, as can be verified from fig 0(b). In this case, the G&P 
formula does not agree with the quantum result. In the computation of the tunneling term 
we have not considered multi-period trajectories, since we have found that the corresponding 
contributions are vanishingly small. 

V. SPIN SYSTEMS 

As discussed in the beginning of section II. c, the coherent state propagator can be de- 
scribed by a variety of semiclassical expressions. For simplicity we chose the one involving 
the classical Hamiltonian 7i. For spin systems there is no direct classical Hamiltonian and 
we have to use a different representation. Our starting point is the spin coherent state, 
defined by \w) = (l + w*w)~ s exp{wS + /H}\s, —s), where Hs is the total spin and w is a com- 
plex number. The Q-representation of the Hamiltonian operator is defined as 7i = (w\7i\w). 
This function appears naturally in the first form of path integral suggested by Klauder and 
Skagerstan [12], whose semiclassical limit requires the addition of an extra term in (|T3|) jl^ . 
the so-called Solari-Kochetov phase [3, 2oT ]: 

"* r d f(l + w*w) 2 dH\ d f(l + w*w) 2 dU 



o 



+ 



At' . (58) 



dw* \ 2s dw J dw \ 2s dw* 
Furthermore, since w and w* are not canonically conjugated, it is convenient to define new 



canoncial variables by 



21] 



Q = H T— ; (w*-w), P= -——( w * + w ). (59) 
V 1 + w*w V 1 + w w 

This transformation is the classical analogue of the Holtein-Primakoff bosonization procedure 



22j. In these variables, the equations satisfied by the stationary phase trajectories that 
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contribute to the semiclassical propagator are the usual Hamilton's equation of motion 

dQ OH dP &H 



(60) 



dt dP ' dt dQ 
We remark that the Solari-Kochetov term does not change the stationary exponent con- 
dition since it already involves second derivatives of the Hamiltonian. With this procedure 
the formalism we have developed can be readily applied to spin systems. The semiclassical 
partition function reads 

i rati) I T (i) \ 1 (-eU)jJfU) \ 

z sc = I . = > or z sc = I . = ' ( 61 ) 

i yTr[M^]-2 j yTr[M 0) ]-2 

for transformations (JSJ and (j2J, respectively. 

As a simple example we consider a spin in a uniform magnetic field. The Hamiltonian 
operator is H = u)S z , leading to the exact partition function 

2s „fej/3(s+l/2) _ -^(8+1/2) 

Z = e^f Y e-^ n = „ fl/2 - . (62) 

n=0 

From the definition of spin coherent states we obtain 

H = -hws =fa(^- + ^-nu,a t (63) 

from which one sees that the only contributing orbit is the trivial one. This yields iS/H = 
Hujfls and iTsx/h = fuuif3/2. The corresponding semiclassical partition function is 

e nu,f3(s+l/2) 
^ sc = e fiL0f3/2 _ e -Kwf3/2 ' 

which differs from the exact result by a factor that is exponentially small in the low temper- 
ature and large spin limits. At first glance it might be expected that the semiclassical result 
should be exact, since the effective Hamiltonian is quadratic. However, it is important to 
note that the classical problem we are dealing with is not completely equivalent to that of 
a harmonic oscillator. In the present case, the phase-space is compact since 7i < fiuis ( see 
Eq. ([63)1 ). The two problems are strictly equivalent only if s — > oo, when the semiclassical 
result is exact. 



VI. CONCLUDING REMARKS 



We have obtained a semiclassical trace formula for the canonical partition function in the 
low temperature limit that does not require sampling and numerical integration in phase- 
space. The examples presented in sections III and IV show that, in spite of the failure of the 
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classical partition function in describing thermal effects at low temperatures, the appropriate 
combination of classical ingredients appearing in the semiclassical formula is indeed capable 
of giving very accurate results in this limit. The quantum corrections come precisely in the 
way we combine the classical functions. The two and three-dimensional extensions of the 
formalism are straightforward in the case of single minimum potentials. 

In the case of multiple minima, the calculation of the non-trivial periodic orbits is more 
involved due to the higher dimensionality. However, contrary to the situation in usual trace 
formulas, only very long orbits should contribute in the low temperature regime. In particu- 
lar, heteroclinic orbits connecting the top of the inverted wells should be of great importance. 
Calculations for two-dimensional potentials are currently under way. Other potentially in- 
teresting perspectives are also open, such as the inclusion of spin-orbit interactions, whose 
semiclassical propagator has been recently derived 23^ . 
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APPENDIX A: EXPANSION COEFFICIENTS 

In this appendix we indicate how to obtain the equalities (|55|). The first expression is 
immediate 

*=VE =«■-*• (A1) 

where we have used the second relation in (J3HJ). Now let us calculate A( u y We have 

From equation (}4*2*|) we have M vv = (1 + 2ujtu'v ') e Qr , which in terms of the independent 
variables v! and v" becomes M vv = (e Qr + 2utu'v"). Therefore 

~M= 2UTV +T d^ 6 ■ (A3) 

One can implicitly determine the derivative of Q starting from Q = 2u(u'v' + a) = 
2u(u'v"e- nr + a). The result is 

~du' ~ 1 + 2utu'v" e- Qr ^ \£m) ~ ' ^ ' 
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Consequently we have ( a ^, v ) = and A( u ) = 0. A similar calculation shows that = 
For A( uv ) we get 

1/ 1 d 2 M vv Z_d^dM^\ 1 / 1 d 2 M vv \ 

{UV) " 2\ M U 2 9u'dv" 2MU 2 du> dv " A 2 \M^ 2 du'dv") Q ' 1 ' 

with 



^)r r (^' e ^ +2 ^ (A6) 



From equation (|A4J) we can show that 



- (Ski-**. (AT) 



du'dv" J \ du'dv" 

which leads to = — 2u;r e _3war . Using the same strategy one can prove that A/ uu ) = 

^.(ot) — 0, that all third order derivatives of $ vanish at the trivial orbit, and that the only 
non- vanishing fourth order derivative is 

*<-KIS)o =4iSwTe ~ w (A8) 
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FIG. 1: Scaled Helmholtz free energy (a) and scaled internal energy (b): higher order semiclassical 
(black), quantum (gray), classical (crosses), and quadratic semiclassical (dashed line). 
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FIG. 2: Scaled entropy: higher order semiclassical (black), quantum (gray), classical (crosses), 
and quadratic semiclassical (dashed line). The third law of thermodynamics is satisfied by the 
semiclassical curves. 
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q T (a.u.) 

FIG. 3: The deep double- well is shown together with the first eleven energy levels in (a). In (b) 
the exact, semiclassical, and G&P results are compared. 
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FIG. 4: The shallow double- well is shown together with the first eleven energy levels in (a). In (b) 
the exact, total semiclassical, harmonic contribution, tunneling contribution, and G&P results are 
displayed. 
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